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PREFACE 

This  Memorandum,  prepared  as  part  of  RAND's  continuing  study  of 
long-distance  optical  propagation  through  the  atmosphere,  gives  a  work¬ 
able  solution  to  the  second-order  Rytov  approximation. 

Recently  there  has  been  a  controversy  ~>ii  the  ranges  for  which  the 
first-order  Rytov  approximation  gives  an  adequate  engineering  estimate 
of  tho  effect  of  a  turbulent  medium  on  laser  propagation.  Therefore, 
it  is  felt  that  an  explicit  derivation  of  the  second-order  Rytov  approx¬ 
imation  should  be  of  interest  to  those  concerned  with  optical  communi¬ 
cation  and  laser  radar.  Furthermore,  the  second-order  approximation 
is  necessary  because  it  is  the  lowest  order  nontrivial  approximation 
which  conserves  average  energy  to  the  order  of  the  approximation. 
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SUMMARY 

\  t 

An  explicit  and  useful  formulation  of  the  solution  for  the 
second-order  Rytov  approximation  is  given.  From  this  solution  a  con¬ 
dition  of  validity  for  the  Rytov  solution  is  obtained.  It  is  concluded 
that,  in  general,  both  the  Born  and  Rytov  approximations  have  ie  same 
domain  validity, 
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I.  INTRODUCTION 

A  great  deal  of  theoretical  effort  has  been  given  to  calculations 

(1-3) 

based  on  the  so-called  Rytov  approximation.  However,  the  valid¬ 

ity  of  the  Rytov  approximation  has  been  questioned  by  a  number  of 
authors.  In  a  recent  paper  Fried ^  continues  to  take  issue  with 

the  arguments  of  Hufnagel  and  Stanley^^  and  Brown,  who  claim  t_nat 
the  Born  and  Rytov  approximations  have  the  same  domain  of  validity. 

In  view  of  this  continuing  controversy,  it  is  felt  that  an  explicit 
derivation  should  be  given  of  the  second-order  Rytov  approximation 
(defined  below).  In  this  Memorandum  the  hierarcny  of  Rytov  equations 
is  presented.  The  explicit  solution  (in  quadrature)  of  the  second-or¬ 
der  Rytov  approximation  is  obtained.  In  Ref.  6  the  Born  approximation 
was  examined  in  detail,  and  it  was  assumed  that  the  domain  of  validitv 
of  the  Rytov  solution  was  the  same.  Here  we  present  an  explicit  demon¬ 
stration  that  this  is  indeed  the  case. 

The  solution  obtained  is  correct  through  terms  of  second  order  in 
n^ ,  where  n^  is  the  fluctuating  part  of  the  index  of  refraction.  From 
this  solution  a  condition  of  validity  is  obtained  for  the  Rytov  approx¬ 
imation  with  the  conclusion  that  both  the  Born  and  Rytov  approximations 
have  the  same  domain  of  validity. 

As  in  Ref.  6,  the  analysis  is  time  independent  and  restricted  tc 
weakly  inhomogeneous  media  which  are  assumed  to  be  statistically  homo¬ 
geneous  and  isotropic,  characterized  by  an  index  of  refraction  corre¬ 
lation  function.  Furthermore,  the  electrical  conductivity  and  magnetic 
permeability  of  the  medium  are  taken  as  zero  and  one,  respectively. 


Scalar  waves  at  optical  wavelengths  arc  considered,  with  the  extensions 


to  vector  fields  being  straightforward.  For  further  details  of  the 
tneory  of  optical  propagation  through  a  turbulent  medium,  the  reader 
is  referred  to  Ref.  6  and  the  refere  ices  therein. 

The  scalar  wave  equation  is 

V2U  +  k2n2(r)U  =  0  (1) 

where  U  is  a  typical  component  cf  the  field,  k  is  the  optical  wave 

number ,  and 


n(r)  =»  1  +  n  .  (r)  (2) 

1 

where  n^  is  the  fluctuating  part  of  the  index  of  refraction  assumed 
to  satisfy  |n^|  s.<  1  and  =  0  (the  bar  over  a  quantity  indicates  the 


ensemble  average  of  the  quantity). 
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II.  THE  SECOND-ORDER  RYTOV  APPROXIMATION 

The  Rytov  transformation  consists  of  setting  U  "  exp[  !<  ]  in 
Eq.  (1)  and  thus  obtaining 


+  (7*) 2  +  k2  [I  +  n1(r)]2  -  0  (3) 

r  1-3) 

In  the  literature  it  is  customary  to  seek  a  solution  for  ^  as  a 

power  series  in  n^.  Let 


if  -  )  (*) 

m 

171*0 

where  '1q  is  zero  order  in  n^,  i}r ^  is  first  order  in  n^,  and  so  forth. 

Then  it  can  be  shown  that  the  i  ‘s  satisfy 

m 


+  (^Q)2  +  k2  *  0 

(5) 

V2^  +  2V^0  •  +  2k2n1(r)  -  0 

(6) 

72!f0  +  27^0  •  Vh2  +  k2n2(r)  +  (V^^2  -  0 

(7) 

and 

m-1 

A 

m 

+  2vl'n  •  Vt  +  \  Vi  •  W  *  o,  m  *  3,4,5... 

U  m  p  m-p 

P=° 

(8) 

Equation 

(5)  gives  the  solution  to  Eq.  (1)  for  the  case  n^  0. 

Thus 

2  2 

if  Uq  =  exp  -Q,  then  Eq.  (5)  is  equivalent  to  v  Uy  +  k  UQ  *  0. 

The  function  T  is  what  is  usually  referred  to  in  the  literature^ 
as  the  Rytov  approximation.  In  this  Memorandum,  ^  is  referred  to  as 
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the  first-order  Rvtov  approximation,  while  ■■  is  referred  to  as  the 
second-order  Rytov  approximation. 

It  can  be  shown  ^  that  within  an  arbitrary  constant  the  general 
solution  of  Eq.  (6)  is  given  by 


’1 


2k^ 

- -  1  dr'  G(r  -  r ')  n  (r')  iMr') 

UQ(r)  1  0 


(9) 


where  the  integration  in  Eq.  (9)  extends  over  the  region  of  space 
where  n^(r)  is  different  from  zero  and 


G(r  -  r')  -  - 


,ik|r-r'j 
L!r-7'|  ~ 


(10) 


is  the  Green's  function  in  the  absence  of  turbulence,  which  satisfies 
the  outgoing  radiation  condition  at  infinity. 

Others  have  failed  to  obtain  a  useful  workable  solution  to 

(3) 

(7)  (e.g.,  Schmeltzer  ),  It  is  noted  that  Eq.  (7)  can  be  solved 
by  the  use  of  a  transformation  suggested  by  a  comparison  of  the  Born 
and  Rytov  expansions.  Upon  equating  second-order  terms  in  the  expansions 
we  have 

2 

’  2  *  ”  2  (11) 


wh^re  .  „  is  given  in  Eq.  (9).  Upon  substituting  Eq.  (11)  into  Eq. 
(’’)  and  using  Eq.  (6),  we  find 


2  7 


0 


2 


,.2  2 

(k  n^ 


2k2uri) 


(12) 


i 


i 

i 


2<r;  - 


Vr> 


id  substitute  this  into  Eq.  (12).  We  find  that  W  satisfies 


99  2  2  2 

7  W  +  k“W  =  -  (k  +  2k  n^V^)  Uq 


Hence , 


W(r)  =  -  k2  [  G(T  -  r')  j^njfr ')  +  2^(70  i^Oj  UQ(7')  dr'  (14) 


Thus,  the  second-order  Rytov  approximation  is  given  by  Eq.  (11), 
where  t,  and  W  are  given  by  Eqs.  (13)  and  (14),  respectively. 

The  transformation  employed  to  solve  Eq.  (7)  suggests  that  the  solu¬ 
tion  to  Eq.  (14)  should  be  equal  to  the  second-order  Born  approximation. 

't 

The  Born  expansion  for  U(r)  is  given  by 


where 


U(r)  =  )  U  (r) 


'  2  L  ,  2 

V  +  k 


K-° 


Uj(r) 


-  2k  '  1  dr'  G(r  -  r ')  n^r')  UQ(r') 


and  for  i 
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Wr) 


2k'L  j  dr'  G(r  -  r ')  n^r ')  U.(r') 


-k2  r  dr'  G(r  -  r')  nj(r')  L’i_1  Cr  0 


^  IT 


(18) 


For  i  *  1  Eq.  (18)  gives 


U2(r)  - 


~  2k2  f  dr'  G(r  -  r ')  n^r  ')  U^r') 


-  k2  f  dr'  G(r  -  r')  nj(r')  U  (r') 


(19) 


Comparing  Eq.  (14)  with  Eq.  (19),  it  is  seen  that  to  within  an  arbitrary 
constant  W(r)  =  U2(r).  Hence,  to  within  an  arbitrary  constant  the 
second-order  Rytov  approximation  is  given  by 


U2(r)  't-(r)  U2(r) 

*•»<*)  *  — - - — =  ~  -  %— r! 


(20) 


uQ(r) 


«n<0 


'uo(r)J 


L_ 


N* 


The  average  field  is  determined,  to  terms  of  second  order  in  n^ , 


as 


U  -  exp  |  VQ  +  tj  +  v2 


exp 


r*.  ♦  r2i 


(21) 


where  (r)  =  j  U„(r)/U„(r) 1.  For  example,  when  the  field  in  the 
absence  of  turbulence  is  a  plane  wave 


B  (0)  /t2. 


ikR 


.  _  (  nv  ,  /k  \  p  e  ’  '  _  ik*Rl 

ik  •  r  (— j—  +  |  «  —  Bn(R)  e  \ 


(22) 


I  I 
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where  B^CR)  n^r)  (r  +  R)  i3  the  correlation  function  of  the 
index  of  refraction  1 luctuations. 

It  is  shown  in  Ref.  6  that  one  must  include  the  second-order 
contribution  to  the  field  (t'2)  in  order  to  obtain  an  adequate  engi¬ 
neering  approximation  to  the  scattered  field  in  the  presence  of 
turbulence.  Furthermore,  a  second-order  expansion  for  the  field 
is  the  lowest  order  nontrivial  approximation  which  conserves  average 
energy  to  the  order  of  the  «... proximation.  That  is,  the  field 

U(T)  =  exp  M'Q(r)  +  ^(’r)  +  ^(r)]  (23) 

which  contains  all  terms  through  second  order  in  n^  gives  the  correct 
average  field,  correct  phase  and  intensity  statistics,  conserves  aver¬ 
age  energy,  and  satisfies  the  optical  theorem  to  this  order.  In  the 
present  notation  the  optical  theorem  is,  through  terms  of  second 
order  in  n^,  given  by 

2Recc2  +  l^l2  -  0  ( 24)1* 


ic 

In  the  usual  notation  the  optical  theorem  is 
where  c  is  the  total  scattering  cross  section  per 
i(0)  is  the  forward  scattering  amplitude  per  unit 


*  c  -  -jT  Imf(O), 

unit  volume  and 
vol ume . 
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III.  VALIDITY  CONDITION 


To  obtain  the  condition  of  validity  for  the  Rytov  solution 


U  “  exp  r Tq  +  t ^  j  it  is  sufficient  to  require  that  each 


successive  term  in  the  expansion  of  V  be  smaller  than  the  preceding 


term.  In  view  of  this,  and  since  the  various  terms  in  the  expansion 


of  ¥  are  random  functions  which  can  be  described  only  in  statistical 


terms,  it  is  assumed  that  the  condition  of  validity  on  the  Rytov  solu¬ 


tion  is  related  to  the  relative  magnitude  of  the  square  value  of  suc¬ 


cessive  terms.  Therefore,  the  condition  of  validity  for  the  Rytov 


solutions  is  taken  to  be 


o2)2  «  (^)2 


f  1 

where  ^  *  y2  "  “*  Now>  ic  can  be  shown  that  |lm  iRe^j 

for  all  cases  of  interest:  hence,  by  expressing  Re'. ^  In  terms  of  V ^ 


by  using  Eq.  (24),  we  find  from  Eq.  (2b,  that  two  conditions  must 


be  satisfied  (from  tb«  real  and  imaginary  parts,  respectively). 


The  real  part  yields  the  condition  (ReV^)  1>  while  the 


imaginary  part  yields  the  stronger  condition,  since  it  can  be  shown 


that  (ReVj)2'  £  (Im 


jihi2  - 1 


Taking  the  ensemble  average  of  Eq.  (26),  we  conclude  that  the 


condition  of  validity  of  the  Rytov  solution  is  given  by 


See  Ref.  1,  Chap,  7  and  9. 


I*-"*:.  >' 
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It  is  noted  that  when  the  successive  terms  in  the  asymptotic 

expansion  for  '{(  (i.e.,  i'  =  +  i'  +  +  •••)  will  all  be  of  the  same 

order  of  magnitude.  The  solution  given  by  Eq.  (23)  is  therefore  valid 
only  when  Eq.  (27)  is  satisfied. 

The  condition  of  validity  of  the  Born  approximation  is  given  by 
exactly  the  same  expression  (Eq.  (57)  of  Ref.  6). 

We  conclude  that,  in  general,  both  the  Bom  and  Rytov  approxi¬ 
mations  have  the  same  domain  of  validity.  A  critical  range  Rc  is 
determined  by  solving  Eq.  (27)  with  an  equality  sign.  For  R  <  R^, 
the  Born  and  Rytov  solutions  are  valid,  while  for  R  ~  R^  both 
solutions  fail  to  yield  a  valid  solution.  Treatments  which  apply 

the  Rytov  approximation  beyond  R^_ ,  and  therefore  predict  large  dis- 
,  .  •  (7) 

persion  of  intensity,  are  incorrect. 
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